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Abstract 

We have searched for decays of the r lepton into seven or more charged particles, using data collected with the 
OPAL detector from 1990 to 1995 in e'^e" collisions at ^/s « Mz- No candidate events were found and an 
upper limit on the branching ratio for t decays into seven charged particles of 1.8 x 10~^ at the 95% confidence 
level was determined. 
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1 Introduction 



During the years 1990 to 1995 of LEP operation, an integrated luminosity of 163 pb~^ has been recorded with 
the OPAL detector at ^/s « 91 GeV, corresponding to approximately five million observed Z° decays. This 
large data sample, which contains more than 210000 r pair events, allows searches for rare decays of the r lepton 
with branching ratios down to about 10~^. The high- multiplicity decays of the r are particularly interesting 
for a measurement of the r neutrino mass. Recent examples of such measurements using r decays into three 
or five charged particles can be found in references ||^, |^. Taking only kinematical considerations into account, 
the decay of the t into twelve pions is possible. 

In this note we present a search for r decays into seven charged particles (7-prong). The HRS Collaboration 
has published an upper limit of 1.9 x 10~^ on the branching ratio for r decays into seven charged particles]^, 
BR (t~ 4tt~ 37r+ Vt (n > 0)), at the 90% confidence level The CLEO Collaboration has determined 
a 95% confidence level upper limit of 1.1 x lO^"' on the branching ratio BR (r^ — > 3tt~ 27r+ 27r" Vr) 0. 



2 The OPAL detector 



A detailed description of the OPAL detector can be found elsewhere [||. Sub-detectors which are particularly 
relevant to the present analysis are described here briefiy. The central detector consists of a system of tracking 
chambers providing charged particle tracking over 96% of the full solid angle inside a 0.435 T uniform magnetic 
field parallel to the beam axis^. Starting with the innermost components, it consists of a high precision silicon 
microvertex detector which was available from 1991, a precision vertex drift chamber, a large volume jet chamber 
with 159 layers of axial anode wires and a set of z chambers measuring the track coordinates along the beam 
direction. The efficiency for separating hits from two different particles in the jet chamber is approximately 
80% for distances between two hits of 2.5 mm Q and drops rapidly for smaller hit distances. The jet chamber 
also provides energy loss measurements which are used for particle identification. 

A lead-glass electromagnetic calorimeter (ECAL) located outside the magnet coil covers the full azimuthal 
range with excellent hermeticity in the polar angle range of | cos 6*1 < 0.82 for the barrel region and 0.81 < 
I cos^l < 0.984 for the endcap region. The forward detectors and sihcon tungsten calorimeters located at both 
sides of the interaction point measure the luminosity and complete the geometrical acceptance down to 24 mrad 
in polar angle. 



3 Event simulation 



The e+e" t^t^ Monte Carlo events were generated using KORALZ 4.0 0. The dynamics of r decays to 
final states with five or fewer charged particles were simulated with the TAUOLA 2.4 decay library The 
7-prong signal events were simulated by modifying the TAUOLA routine describing the decay of the r into 
five or more pions to allow up to seven charged pions with either zero (r^ Att~ 37r+ Vr) or one neutral pions 
(t~ — > Att~ Stt"*" tt^ Vt) in the final state. To determine the efficiency for reconstructing an event with a 7-prong 
T decay we generated r pair events in which the first r decays according to the standard TAUOLA library and 
the second one decays into seven charged particles. 

We have implemented various phase space models for these 7-prong decays. The standard case assumes 
an isotropic phase space between the minimum allowed value, equivalent to the sum of the pion masses, and 
the maximum allowed value, rriT. To study the effect of the phase space modelling on the final result, we also 
considered the two extreme cases of a "low" and "high" phase space, that is with a 7-pion invariant mass always 

^References in this paper to specific charge states apply to the charge conjugate states as well. 

■^The OPAL coordinate system is a right-handed coordinate system which is defined so that the 2:-axis is in the direction of the 
electron beam, the x-axis is horizontal and points towards the centre of the LEP ring; 9 and rp are the polar and azimuthal angles, 
defined relative to the +z- and -|-x-axes, respectively. 
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below or above the value of 1.4 GeV/c^, which is the mean value of TttItt and m^. Additionally, we generated 
7-prong Monte Carlo events where the invariant mass of the 7-pion system is peaked close to this mean value 
in a Gaussian manner with a standard deviation of 150 MeV/c'^. 

To estimate the background from e^e^ —> t^t^ events with r decays into up to five charged particles, we 
have used a sample of 920000 events generated with KORALZ and TAUOLA. The expected rate of multihadronic 
background was determined using e+e" — > qq events generated with JETSET The background coming 
from e^e~ —^ e+e~ events was estimated using a sample of simulated e"'"e~ — ^ e+e^ events corresponding 
to a luminosity of 287 pb~^ and generated with BABAMC |jl^. Contributions from two-photon processes 
were studied using the PYTHIA |^ and Vermaseren Monte Carlo generators. Background contributions 
coming from e+e~ — > n~ and four-fermion events are small. They were investigated using the KORALZ and 
FERMISV |l^ event generators, respectively. 

All Monte Carlo events were passed through the GEANT [|3) simulation of the OPAL detector Q and 
were then processed in the same way as the real data. 



4 Event preselection 

A cone jet algorithm ||l^ was employed to assign all tracks and electromagnetic clusters to cones with a half 
opening angle of 35°. We required that exactly two cones with at least one track each were found in the event. 
Tracks were required to satisfy the following conditions: 

• p±^ > 0.1 GcV/c, where p±^ is the momentum component transverse to the beam direction. 

• A minimum number of 20 hits in the central jet chamber was required. This restricts the acceptance of 
the detector for tracks to | cos 6*1 < 0.965, where 9 is the angle between the track and the beam axis. 

• The distance of closest approach of the track to the beam axis must be less than 2 cm and the displacement 
of the track along the beam axis from the nominal interaction point at the point of closest approach to 
the beam must be less than 20 cm. 

To reject events coming from two-photon interactions, the acollinearity of the two cones was required to be 
less than 15°, with the direction of each cone given by the momentum sum of all tracks and electromagnetic 
clusters inside the cone. The momenta of the ECAL clusters were calculated using the energy deposition and the 
angles of the cluster. The visible energy of each jet is taken as the maximum of the sum of the track momenta 
and the sum of the ECAL cluster energies associated with that jet. Events were rejected if the sum of the 
visible energies of the two jets was less than 3% of the centre-of-mass energy. Events with a total visible energy 
smaller than 20% of the centre-of-mass energy were rejected if the sum of the transverse momenta of the tracks 
and the sum of the transverse momenta of the electromagnetic clusters were both smaller than 2 GeV. This 
set of cuts rejects most of the two-photon background. Cosmic ray events were removed using the time-of-flight 
and the tracking chamber information ||l6|. Finally, we rejected events where some particles may be lost close 
to the beam pipe by requiring that the total energy deposit in the forward detectors be less than 2 GeV. We 
also required | cos^concl < 0.9 for each of the two cones in the event, where ^cono is the angle between the jet 
direction and the beam axis. This cut ensures that both jets are located inside the acceptance of the detector. 

Only events with a minimum of four tracks were accepted for further analysis. The efficiency for selecting a 
T+r^ event with a 7-prong t decay after the preselection cuts was (80.8 ± 0.5)% according to the Monte Carlo 
simulation. 



5 Selection of 7-prong r events 

At this stage, the largest sources of background, around 55%, come from e+e" — > qq events with a low- 
multiplicity jet, followed by e"'"e~ — > t^t~ events with r decays into five or fewer charged particles (« 33%). 
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These backgrounds were reduced by the following set of cuts. Since approximately 99.9% of the r leptons decay 
into one or three charged particles, we required that the number of tracks iVi in one of the two cones^ is exactly 
1 or 3. The absolute value of the sum of the charges of these tracks was required to satisfy | X^tracks^li = ^■ 
The search for the 7-prong r decay was performed in the second cone. We required that this contains at least 
six tracks {N2 > 6) and that the absolute value of the sum of the track charges satisfies | X^tracks 9 1 2 ^ 2. The 
distribution of N2 is shown in Figure |^. 

The invariant mass of the tracks from a true 7-prong t decay must be less than the r mass. Therefore the 
invariant mass of the tracks in the second cone, m2, calculated using the pion mass for each track, was required 
to satisfy 1112 < 1.6 GeV/c^ for cones with six reconstructed tracks and 777,2 < 1-9 GeV/c^ for cones with 
N2 > 7. The cut is lower for the N2 = 6 case since the measured invariant mass of a true 7-prong decay 
candidate would be smaller due to the missing track. The distribution of 7772 before the cut is shown in Figure ^. 

Most of the qq events that were not removed by the invariant mass cut were rejected by a cut on the 
maximum opening angle of any pair of tracks in the two jets. For the first cone we required this angle to be 
less than 10° (a™'''' < 10°) when A^i = 3. For the second cone we required a™'"'' < 8°. The distribution of 
Q,max jg g]^Q.^^jj jjj Figure || for the data and for the signal and background Monte Carlo samples after the cuts 
on 7772 and af^'^^ have been applied. 

Jets likely to contain one or more electrons coming from r decays with five or fewer tracks and additional 
tracks from photon conversions or 7r° Dalitz decays were rejected using the energy loss measurement and 
information on the energy deposit in the electromagnetic calorimeter. For the tracks in the 7-prong candidate 
cone we have calculated a quantity /2, which is the logarithm of the ratio of the probability for the measured 
dE/dx of the track to be consistent with an electron to the probability for the measured dE/dx of the track to 
be consistent with a pion. If at least 20 dE/dx measurements along the specified track are available, we required 
/2 to be less than 5. For the N2 = Q case this cut was lowered to /2 < 3 to reject the conversion background 
more effectively. No cut was made for tracks with less than 20 dE/dx measurements or for tracks where both 
probabilities were smaller than 1%. The distributions of the maximum /2 is shown in Figure ^. In this figure, 
the expected signal contribution is normalised assuming a branching ratio equal to the upper limit quoted in 
reference § . To reduce further the background from events containing electrons, we applied the additional cut 
(-E'cone/-Pcone)2 < 0.8, whcrc -Econc IS the sum of the energies of all electromagnetic clusters inside the cone and 
fcono is the sum of the momenta of all tracks associated with the cone. 

After all cuts, no event is found in the data. The expected number of background events from all modelled 
sources is O.SlJ'g. The numbers of expected and observed events after each cut are shown in Table 0. The 
efficiencies for selecting a 7-prong t event arc shown in Table ^ for the signal Monte Carlo samples described in 
section ^. A branching ratio equal to the upper limit quoted in reference would result in approximately 33 
T pair events containing a 7-prong t decay being observed in our data sample. 



6 Systematic errors 

The two sources of systematic errors in this analysis are the uncertainty in the total number of r pair events, 
Nt+t- , and the uncertainty in the 7-prong selection efficiency. We determined N^+t- , using the total number of 
hadronic Z'^ decays and dividing by the ratio of the hadronic and the leptonic branching ratio of the Z*' corrected 
for contributions from s-channel photon exchange. The calculated number of t pair events is N-^+r- = 210680± 
1100. This corresponds to a systematic error of 0.5%. 

For the determination of the upper limit on the branching ratio we chose the 7-prong Monte Carlo simulation 
with the isotropic phase space distribution of the 7-pion system. After all cuts we obtain a selection efficiency 
of 0.410 for this 7-prong Monte Carlo sample. Among all other 7-prong Monte Carlo simulations described in 
section only the high phase space model leads to a lower final selection efficiency, which is still consistent with 
the selection efficiency for the isotropic phase space case. As an estimate of the systematic uncertainty from 
phase space effects, we therefore use the statistical error of the selection efficiency obtained for the isotropic 
7-prong Monte Carlo sample, corresponding to a relative error of 3.9%. 

•^In the rest of the paper we refer to this as the first cone. 
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Cut 


Number of expected events from MC 
r+T^ e+e" qq Sum of MC 


OPAL 

data 


Preselection, 
TVi = 1 or 3, 
N2>6 


788±13 


83±7 


8684±79 


9555±81 


9593 


1 Etracks'jll = 1 
lEtracks^b < 2 


657±12 


79±7 


8058±76 


8794±78 


8888 


invariant 
mass cuts 


84.6±4.4 


20.5±3.4 


35.4±5.1 


140.5±7.5 


131 


Q,max ^ go 


50.0±3.4 


18.2±3.2 


0.6±0.6 


68.8±4.7 


67 


dE/dx cuts 


3.0±0.8 


1.7±1.0 


< 1.1 


4 7+1-7 

^■'-1.3 


4 


(-^conc/-fconc ) 2 

requirement 


0.5±0.3 


< 1.1 


< 0.6 








Table 1: Number of expected and observed events for a luminosity of 163 pb~ after application of each cut. 
Only the most important sources of background are included. The number of expected background events is 
obtained using the Monte Carlo samples (MC) described in the text. The errors shown are statistical only. 

To take into account that additional neutral pions may be produced in the final state, we checked differences 
in the selection efficiencies for the 47r~ Stt^ Vr and the 47r~ 37r^ tt*^ Vr decays. Since we did not obtain a lower 
selection efficiency for the 47r~ 37r''" tt" Vr case, we conclude that the branching ratio limit is relatively insensitive 
to the number of additionally produced neutral pions. The case that kaons may be produced in the final state 
has been neglected because the corresponding decay modes are expected to be highly suppressed compared to 
the modes we have considered. 

We studied the possible differences in the modelling of the number of reconstructed tracks between data and 
the Monte Carlo simulation. This comparison was performed using the dE/dx distribution for the tracks in 
the 7-prong candidate cone. In the OPAL jet chamber, for momenta above 2 GeV/c, single tracks should have 
a maximum measured d£'/da; of 10 keV/cm (electrons) and only unresolved double tracks should have a larger 
dE/dx value. If at least one track with a momentum greater than 2 GeV/c and a dE/dx value greater than 
15 keV/cm was found, we conclude that this is due to two tracks that have not been resolved. After applying 
the preselection and the multiplicity cuts we calculated the fraction /ics of 7-prong candidate cones in which 
such a high dE/dx value was measured. We found that f^^s is about 2% smaller in the data than in the Monte 
Carlo simulation. This 2% discrepancy is also found in the low invariant mass region selected with TO2 < 2 GeV. 
Relaxing the cut on iV2, we repeated this study for A'2 = 4 and N2 = 5 and obtained similar results. Because 
a higher efficiency for the data than for the Monte Carlo simulation leads to a systematic error which does not 
affect the final limit, this effect has been neglected. 

The uncertainty arising from the cut on the total charge of the cones was estimated by comparing the ratio 
of rejected events for data and Monte Carlo simulation when the cut is applied directly after the preselection 
and the multiplicity cuts, giving a relative uncertainty of 0.7%. For the invariant mass cut and the cut on the 
maximum angle between any pair of tracks, the systematic uncertainties were estimated by comparing data and 
Monte Carlo samples selected by requiring after the preselection iVi = 1 and 5 < iV2 < 7. The distribution of 
TO2 after applying these selection criteria is shown in Figure ^. For the two quantities, we obtain a combined 
relative systematic error of 1.1%. To determine the systematic uncertainties introduced by the cuts on /2 and 
(-E'conc/-Pconc)2, wc havc Compared data and Monte Carlo samples selected by requiring after the preselection 
cuts Ni — 1 and N2 > 6 with an additional cut on the invariant mass in the second cone of TO2 < 3 GeV. For 
these samples, containing mostly r pair events, a comparison of the /2 and (i?conc/^conc)2 distributions between 
data and the Monte Carlo simulation leads to 7-prong efficiency uncertainties of 0.6% and 1%, respectively. 
These results are also consistent with those obtained in ||l^ . 

Adding the systematic errors in quadrature, the total relative uncertainty of the 7-prong selection efficiency 
is 4.3% and is dominated by the uncertainty in the modelling of the 7-prong phase space. 
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Cut 


Efficiencies for the 7-prong selection 


At: 37r+ Vr 

isotropic peak high low 


An 37r+ t:^ Vr 
isotropic 


Preselection, 
TVi = 1 or 3, 
TVs > 6 


0.612±0.015 


0.630±0.015 


0.660±0.015 


0.627±0.015 


0.645±0.011 


1 Etracks^ll = 1 
1 Etracks^b < 2 


U.O / ZztU.UiO 


U.DU4ztU.UiO 


U.OloztU.UiO 


u.oyzztu.uiD 


U.DUoztU.Ull 


invariant 
mass cuts 


0.472±0.016 


0.503±0.016 


0.453±0.016 


0.513±0.016 


U . tJ Z O ZC U . U J- z 


a^'^'^ < 8° 


0.447±0.016 


0.480±0.016 


0.444±0.016 


0.479±0.016 


0.493±0.012 


dE/dx cuts 


0.419±0.016 


0.453±0.016 


0.416±0.016 


0.451±0.016 


0.443±0.012 


(-^cone / -Pcone ) 2 

requirement 


0.410±0.016 


0.446±0.016 


0.404±0.016 


0.445±0.016 


0.420±0.012 



Table 2: Efficiencies after each cut for the selection of a 7-prong t event as obtained from the different Monte 
Carlo simulations described in section The errors shown are statistical only. 



The reliability of the predicted background rate for e+e" e+e" events and c+c~ t'^t~ events with t 
decays into up to five charged particles was effectively tested using the distribution of the total electromagnetic 
energy deposited in the second cone before the dE/dx cuts, as shown in Figure^. Good agreement is observed 
between the data and the Monte Carlo simulation. The number of expected multihadronic background events 
has been obtained using the JETSET prediction which has been scaled down by 10-20% (depending on the 
exact fragmentation parameters and hadronic branching ratios used in each JETSET sample) to fit the data for 
events with iVi = 1 or 3 and 10 < N2 < 15. To check the reliability of the e+e" qq background prediction, we 
have repeated the analysis without applying the dE/dx and (-Econo/Jconc)2 cuts and reversing the cut on a^'^^. 
This selects an almost pure sample of multihadronic decays in the region where the signal is expected. Within 
the limited statistics available, good agreement is found between the data and the Monte Carlo prediction: 5 
events are observed in data and 4.8 ±1.8 are predicted by the simulation. 



7 Results 



Using data collected with the OPAL detector at LEP at centre-of-mass energies on or near the Z° resonance, 
we have searched for decays of the r lepton into seven or more charged particles. No candidate events were 
observed. An upper limit on the branching ratio is obtained using the calculated number of r pair events, 
N^+r- = 210680 ± 1100, and the detection efficiency of (41.0 ± 1.8)% according to 

3 

BR (t- Att- 3n+ n'j Vr (n > 0)) < '■ . (1) 

^ / r V - 2-iV^+^- -0.410 ^ ' 

The total relative systematic error of 4.3% was included in the manner described in reference iQ. We obtain 
the upper limit on the branching ratio of 

BR (j- An- 3n+ ^7 j/^ (n > 0)) < 1.8 x lO"'^ (2) 

at the 95% confidence level. This is equivalent to an upper limit of 1.4 x 10^^ at the 90% confidence level. 
An improvement of more than one order of magnitude is achieved compared to the only previously published 
result 01. 
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Figure 1: Distribution of the number of tracks, N2, in the second cone after relaxing the cut to N2 > 5. An 
arbitrary normalisation is used for the 7-prong signal Monte Carlo simulation. The other histograms show the 
contributions from the various background Monte Carlo samples. The points are the data. The arrow indicates 
the position of the selection cut. 
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Figure 2: Distribution of m2 before the invariant mass cut (a) and of a^^^ after the cuts on m2 and a™^^ 
(b). An arbitrary normalisation is used for the 7-prong signal Monte Carlo simulation. The other histogram,s 
show the contributions from the various background Monte Carlo samples. The points are the data. The arrows 
indicate the positions of the selection cuts. 
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Figure 3: Distribution of the maximum /2 before the dE/dx cuts. The 7-prong signal Monte Carlo simulation 
was normalised assuming a branching ratio equal to the upper limit quoted in Ref. [4]- The other histograms 
show the contributions from the various background Monte Carlo samples. The points are the data. The arrows 
indicate the position of the selection cuts. 
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Figure 4: Distribution of m2 selected by requiring after the preselection Ni ~ 1 and 5 < N2 < 7 ('aj 
and of the total electromagnetic energy deposited in the 7-prong candidate cone before the dE/dx cuts (h). An 
arbitrary normalisation is used for the 7-prong signal Monte Carlo simulation. The other histograms show the 
contributions from the various background Monte Carlo samples. The points are the data. 
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